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'~~ Abstract. The purpose of the present manuscript is to collect known results 

pJT and present some new ones relating to nodal domains on graphs, with special 

^^H emphasize on nodal counts. Several methods for counting nodal domains will 

On be presented, and their relevance as a tool in spectral analysis will be discussed. 
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1. Introduction 



Spectral graph theory deals with the spectrum and the eigenfunctions of the 
Laplace operator defined on graphs. The study of the eigenfunctions, and in partic- 
ular, their nodal domains is an exciting and rapidly developing research direction. 
It is an extension to graphs of the investigations of nodal domains on manifolds, 
CO which started already in the 19th century by the pioneering work of Chladni on the 

nodal structures of vibrating plates. Counting nodal domains started with Sturm's 
oscillation theorem which states that a vibrating string is divided into exactly n 
^b nodal intervals by the zeros of its n*'' vibrational mode. In an attempt to general- 

(^ ize Sturm's theorem to manifolds in more than one dimension, Courant formulated 

' his nodal domains theorem for vibrating membranes, which bounds the number 

^—{ of nodal domains of the n*'' eigenfunction by n [l] . Pleijel has shown later that 

t"~-- Courant 's bound can be realized only for finitely many eigenfunctions [2]. The 

^D study of nodal domains counts was revived after Blum et al have shown that nodal 

^ count statistics can be used as a criterion for quantum chaos [3] . A subsequent pa- 

per by Bogomolny and Schmit illuminated another fascinating connection between 
nodal statistics and percolation theory [4]. A recent paper by Nazarov and Sodin 
5^ addresses the counting of nodal domains of eigenfunctions of the Laplacian on §^ 

[5] , They prove that on average, the number of nodal domains increases linearly 
with n, and the variance about the mean is bounded. At the same time, it was 
shown that the nodal sequence - the sequence of numbers of nodal domains ordered 
by the corresponding spectral parameters - stores geometrical information about 
the domain J6] . Moreover, there is a growing body of numerical and theoretical 
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evidence which shows that the nodal sequence can be used to distinguish between 
isospectral manifolds. O [8]. 

In the present paper we shall focus on the study of nodal domains on graphs, 
and show to what extent it goes hand in hand or complements the corresponding 
results obtained for Laplacians on manifolds. 

The paper is designed as follows: The next chapter summarizes some elemen- 
tary definitions and background material necessary to keep this paper self contained. 
Next, we survey the known results regarding counting nodal domains on graphs and 
state a new theorem regarding the morphology of nodal domains. After these pre- 
liminaries, we present a few counting methods of nodal domains on graphs. Finally, 
the intimate connection between nodal sequences and isospectrality on graphs will 
be reviewed, and some open problems will be formulated. 



2. Definitions, notations and background 

A graph Q = (V,S) is a set of vertices V = {1,2, ■ ■ -V} of size V = |V| and a set 
of undirected bonds (edges) B of size B = \B\, such that {i,j} £ B ii the vertices 
i and j are connected by a bond. In this case we say that vertices i and j are 
adjacent and denote this by ? ~ j. The degree (valency) of a vertex is the number 
of bonds which are connected to it. A graph is called ii-regular if all its vertices 
are of degree v. Throughout the article, and unless otherwise stated, we deal with 
connected graphs with no multiple bonds or loops (a bond which connects a vertex 
to itself). A well known fact in graph theory is that the number of independent 
cycles in a graph, denoted by r is equal to: 

(1) r^ B-V + Co 

where Co is the number of connected components in Q. We note that r is also the 
rank of the fundamental group of the graph. A tree is a graph for which r — 0. 
Let g he a, subgraph of G- We define the interior of g as the set of vertices whose 
adjacent vertices are also in g. The boundary of g is the set of vertices in g which 
are not in its interior. 

A graph G is said to be properly colored if each vertex is colored so that adjacent 
vertices have different colors. G is k-colorable if it can be properly colored using k 
colors. The chromatic number x{Q) is fc if t/ is fc-colorable and not (A;-l)-colorable. 
A very simple observation, which we will use later, is that x{Q) ^ ^• 

G is called bipartite if its chromatic number is 1 or 2. However, since a chro- 
matic number 1 corresponds to a graph with no bonds, and we are dealing only 
with connected graphs, we can exclude this trivial case and say that for a bipartite 
graph, X = 2. The vertex set of a bipartite graph G can be partitioned into two dis- 
joint sets, say Vi and V2, in such a way that every bond of G connects a vertex from 
Vi with a vertex from V2. We then have the following notation: G = (Vi U V2,B) 

m- 

The adjacency (connectivity) matrix of G is the symmetric V x V matrix C = 
C[G) whose entries are given by: 

„ _ J 1, if i and j are adjacent 
*-' 1 0, otherwise 
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Laplacians on graphs can be defined in various ways. The most elementary way 
relies only on the topology (connectivity) of the graph, and the resulting Laplacian 
is an operator on a discrete and finite-dimensional Hilbert space. These operators 
or their generalizations to be introduced below will be referred to as "discrete" or 
"combinatorial" Laplacians. One can construct the Laplacian operator as a differ- 
ential operator if the bonds are endowed with a metric, and appropriate boundary 
conditions are required at the vertices. The resulting operator should be referred 
to as the "metric" Laplacian. However, because the metric Laplacian is identical 
with the free Schrodinger operator (i.e. with no potential) on the graph, one often 
refers to this system as a "Quantum Graph" - a misnomer which is now hard to 
eradicate. In the sequel we shall properly define and discuss the relevant versions 
of Laplacians on graphs. 

The discrete Laplacian, of Q, is the matrix 

(2) LiG) = D-C , 

where D is the diagonal matrix whose i*'* diagonal entry is the degree of the vertex 
i, and C is the adjacency matrix of Q. A generalized Laplacian, L is a symmetric 
V xV matrix with off-diagonal elements defined by: ij < if vertices i and j are 
adjacent, and L^ ■ = otherwise. There are no constraints on the diagonal elements 
of L'. 

The eigenvalues of L(Q) together with their multiplicities, are known as the 
spectrum oi G- To each eigenvalue corresponds (at least one) eigenvector whose 
entries are labeled by the vertices indexes. It is well known that the eigenvalues 
of the combinatorial Laplacian are non-negative. Zero is always an eigenvalue 
and its multiplicity is equal to the number of connected components of Q. An 
important property regarding spectra of large w-regular graphs is that the limiting 
spectral distribution is symmetric about X = v, and is supported on the interval 
[v - 2y/u^^, V + 2\/i7^^] [18] . 

An extensive survey of the spectral theory of discrete Laplacians can be found 
in [2111221123] . 

To define quantum graphs a metric is associated to Q. That is, each bond 
is assigned a positive length: Li, e (0,cx)). The coordinate along the bond b is 
denoted by Xb. The total length of the graph will be denoted by £ = J2beB-^b- 
This enables to define the metric Laplacian (or free Schrodinger operator) on the 

,2 

graph as the negative second derivative —f;^ on each bond. The domain of this 
operator on the graph is the space of functions which connect continuously across 
vertices and which belong to the Sobolev space W^'^{b) on each bond b. Moreover, 
vertex boundary conditions are imposed to render the operator self adjoint. We 
shall consider in this paper the Neumann and Dirichlet boundary conditions: 

(3) Neumann condition on the vertex i : 

(4) Dirichlet condition on the vertex i : 

where S^'"' denotes the group of bonds which emerge from the vertex i and the 
derivatives in (pi) are directed out of the vertex i. The eigenfunctions are the 



^ dXb '^'^'''^ 


-- 


= 


beSi'i ^^^ -->=" 




Mxb%,=o-- 


= 0, 





4 RAM BAND*, IDAN OREN" AND UZY SMILANSKY"'^ 

solutions of the bond Schrodinger equations: 

(5) ybeB -^V6 = fc'V'b, 

which agree on the vertices and satisfy at each vertex boundary conditions of the 
type ^ or Q. The spectrum {k^}^=i is discrete, non-negative and unbounded. 
One can generahze the metric Laplacian by including potential and magnetic flux 
that are defined on the bonds. Other forms of boundary conditions can also be 
used. However, these generalizations will not be addressed here, and the interested 
reader is referred to two recent reviews |13L 114] . 

Finally, two graphs, Q and Ti, are said to be isospectral if they posses the same 
spectrum (same eigenvalues with the same multiplicities). This definition holds 
both for discrete and quantum graphs. 

3. Nodal domains on graphs 

Nodal domains on graphs are defined differently for discrete and metric graphs. 

• Discrete graphs: Let G — (V, B) be a graph and let f=(/i, /2, . . . , fv) be a real 
vector. We associate the real numbers fi to the vertices of G with i — 1,2, ... ,V. 
A nodal domain is a maximally connected subgraph of G such that all vertices have 
the same sign with respect to f . The number of nodal domains with respect to a 
vector f is called a nodal domains count, and will be denoted by i'{i). The maximal 
number of nodal domains which can be achieved by a graph G will be denoted by 
vg. The nodal sequence of a graph is the number of nodal domains of eigenvectors 
of the Laplacian, arranged by increasing eigenvalues. This sequence will be denoted 

byi^^JLi- 

The definition of nodal domains should be sharpened if we allow zero entries in f . 

Two definitions are then natural: 

• A strong positive (negative) nodal domain is a maximally connected sub- 
graph HofG such that fi>0 (/, < 0) for all i G H. 

• A weak positive (negative) nodal domain is a maximally connected sub- 
graph HoiG such that /^ > (/^ < 0) for all i G H. 

In both cases, a positive (negative) nodal domain must consist of at least one 
positive (negative) vertex. According to these definitions, it is clear that the weak 
nodal domains count is always smaller or equal to the strong one. 

• Metric graphs: Nodal domains are connected domains of the metric graph where 
the eigenfunction has a constant sign. The nodal domains of the eigenfunctions are 
of two types. The ones that are confined to a single bond are rather trivial. Their 
length is exactly half a wavelength and their number is on average — . The nodal 
domains which extend over several bonds emanating from a single vertex vary in 
length and their existence is the reason why counting nodal domains on graphs is 
not a trivial task. The number of nodal domains of a certain eigenfunction on a 
general graph can be written as 

(6K = ^E E { ^ +^(l-(-l)L^Jsign[<^.]sign[^,])|-B + V^ 

where [x\ stands for the largest integer which is smaller than x, and (j>i, 4>j are the 
values of the eigenfunction at the vertices connected by the bond b — {i,j} |17j . 
mj holds for the case of an eigenfunction which does not vanish on any vertex: 
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Vj (j)i ^ 0, and there is no cycle of the graph on which the eigenfunction has a 
constant sign. The last requirement is true for high enough eigenvalues where half 
the wavelength is smaller than the length of the shortest bond. This restriction, 
which is important for low eigenvalues, was not stated in |17j . 

Nodal domains on quantum graphs can be also defined and counted in an 
alternative way. Given an eigenfunction, we can associate to it the vector (p = 
{(pi , . . . , 4>v) of its values on the vertices and count the nodal domains of this vector 
as in the case of a discrete graph, explained above. The reasoning behind this way 
of counting is that the values of the eigenfunction on the vertices {0i}lLi together 
with the eigenvalue k^ store the complete information about the values of the 
eigenfunction everywhere on the graph. We thus have two independent ways to 
define and count nodal domains on metric graphs. To distinguish between them we 
shall refer to the first as metric nodal domains, and the number of metric domains 
in the n*^ eigenfunction will be denoted by fin- The domains defined in terms of 
the values of the eigenfunction on the vertices will be referred to as the discrete 
nodal domains. The number of the discrete nodal domains in the n*'' eigenfunction 
will be denoted by Vn, similar to the notation of this count for the discrete graphs. 

As far as counting nodal domains is concerned, trees behave as one dimensional 
manifolds, and the analogue of Sturm's oscillation theory applies for the eigenfunc- 
tions of the discrete [9] and the metric Laplacians |10l IllL 1121 115] , as long as the 
eigenvector (or the eigenfunction) does not vanish at any vertex. Thus we have 
Vn = n for discrete tree graphs and fin — n for metric ones. 

Similarly, Courant's theorem applies for the eigenfunctions of both the discrete 
and the metric versions of the Laplacian on any graph: Vn < "-, A*™ < n, |16|i I17j . 
It should be noted that there is a correction due to multiplicity of the n*'' eigenvalue 
and the upper bound becomes n + rn— 1, where m is the multiplicity |16j . However, 
sharper lower and upper bounds for the number of nodal domains were discovered 
recently. Berkolaiko provided a lower bound for the nodal domains count for both 
the discrete and the metric cases |27j . He showed that the nodal domains count of 
the n"* eigenfunction of the Laplacian (either discrete or metric) has no less than 
n — r nodal domains (r is the number of independent cycles in the graph). Again, 
this is valid if the eigenfunction has no zero entries and it belongs to a simple 
eigenvalue. When n — r < 0, this result is trivial since a nodal domains count is 
positive by definition. We note that for metric graphs this theorem does not hold 
when the discrete count is used. This can be explained by the simple observation 
that n — r grows unbounded while the discrete count is bounded by the number of 
vertices. 

A global upper bound for the nodal domains count of a graph Q was derived 
in [28j : The maximal number of nodal domains on Q was proven to be smaller or 
equal to j/gV — x + 2, where x is the chromatic number of Q. This bound is valid 
for any vector, not only for Laplacian eigenvectors. 

To end this section we shall formulate and prove a few results which show that 
not all possible subgraphs can be nodal domains of eigenvectors of the discrete 
Laplacians of w-regular graphs. The topology and connectivity of nodal domains 
are restricted, and the restrictions depend on whether the eigenvalue is larger or 
smaller than the spectral mid-point v. 



Theorem 3.1. Let Q be a v-regular graph. Then the following statements hold: 
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i. For all eigenvectors with eigenvalue X > v the nodal domains do not have interior 

vertices. 

ii. For all eigenvectors with eigenvalue X < v, all the nodal domains consist of at 

least two vertices. 

Hi. For all eigenvectors with eigenvalue X < v — k (and k < v), in every nodal 

domain there exists at least one vertex with a degree (valency) which is larger than 

k. 

Proof. Let f be an eigenvector with no zero entries of the discrete Laplacian, 
corresponding to an eigenvalue A. Let 5 be a nodal domain of f . 
i. Assume that i is an interior vertex in g. Hence, the signs of fj for all j ^ i are 
the same as the sign of fi. This is not compatible with 

(7) -^/, = (A-«)/, 

for X > V. Hence g cannot have any interior vertices. 

ii. Assume that the subgraph g consists of a single vertex i. Thus on all its neighbors 

j ^ i, the sign of fj is different from the sign of /«. This is not compatible with 

(8) -J2fl = -iv-X).f^ 

for X < V. Hence g cannot consist of a single vertex. 

iii. Denote the complement of the nodal domain (subgraph) (7 in V by g'^. For all 

the vertices i in g 

(9) m, = Vf, - Y, C,^^f, - Y, Cl,^fl = A/, . 

Summing over i Cz g we get: 

(10) (« - ^) E /^ = E E ^^^-/■^ + E ^m/; 

ieg ieg \jeg leg" 

Assuming for convenience that fi are positive for i (z g, the rightmost sum in the 
equation above is non positive, and therefore 

(11) (« ~ A) ^ /. < YT. ^i^^f^ = E "-/■« ^ ^ E /- 

ieg ieg jeg ieg ieg 

Here, Vi — X^ieq ^ji i^ ^^'^ valency (degree) of the i*'' vertex in g, and v denotes 
the largest valency in the subgraph. Since it is assumed that A < w — fc we get 

(12) k < V , 
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which completes the proof. D 

This theorem holds also for an eigenvector f which has zero elements with the 
only exception being the failure of part i when using the weak count. The case 
X = V deserves special attention. As long as the nodal domain under study has no 
vanishing entries, it cannot consist of a single vertex nor can it have interior vertices. 
Namely, for X — v, statements i and ii of Theorem 3.1. are valid simultaneously. 
Otherwise, one should treat separately the strong and the weak counts. For the 
strong count, and A = w a nodal domain cannot have an interior vertex. However, 
using the weak count for X ~ v one finds that no single vertex domains can exists, 
as in Theorem S.l.ii. 

Item in of Theorem 3.1 can be used to provide a A dependent bound on the 
number of nodal domains of eigenvectors corresponding to eigenvalues X < v. De- 
fine the integer k as k = v — [A] . Theorem S.l.iii implies that every nodal domain 
occupies at least k + 2 vertices. Thus, their number is bounded by -j^. Courant 
theorem guarantees that the number of nodal domains is bounded by the spectral 
count A/'(A). This information together with the known expression for the expecta- 
tion value of A/'(A) over the ensemble of random graphs, enable us to show that for 
large v and V, the bound j^ is more restrictive than the Courant bound. Unlike 
Pleijel's result, this bound is not uniform for the entire spectrum, and it applies 
only to the lower half of the eigenvalues with X < v. 

Theorem 3.1 can be easily extended to the nodal properties of the eigenvectors 
of the generalized Laplacian, provided that the weights at each vertex sum up to a 
constant v which is the same for the entire graph. 



4. How to count nodal domains on graphs? 

When discussing nodal domains counting, we must make a clear distinction 
between algorithmic and analytic methods. In the first class, we include computer 
algorithms. They vary in efficiency and reliability, but they have one feature in 
common, namely, that the number of nodal domains is provided not as a result of a 
computation, but rather, it follows from a systematic counting process. The most 
widely used method is the Hoshen-Kopelman algorithm (HK) for counting nodal 
domains on 2-dimensional domains |30j . Analytic methods provide the number 
of nodal domains as a functional of the function and the domain under study. 
The functional might be quite complicated, and not efficient when implemented 
numerically. An example of an analytical method for nodal domains counting in 
one dimension, is given by 

(13) .= tsifix)) ^^(^) 



dx 



dx 



where the nodal domains of f{x), in the interval [a, b] are provided (assuming that 
f{a)f{b) y^ 0). While counting in 1-d is simple, there is no analytic counting method 
for computing the number of nodal domains in higher dimensions: the complicated 
connectivity allowed in high dimensions renders the counting operation too non 
local. 

Graphs, which are in some sense intermediate between one and two dimensions 
still allow several analytic counting methods which we discuss here. An example 
of an analytical count is given by ^ . The HK algorithm is well suited for graphs 
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which are grids. However, it is not as efficient when the graph under study is 
highly irregular. Although the HK algorithm fails for very complex graphs, other 
algorithms, called labeling algorithms, display linear efficiency ([46], [47'). 
Method ///. in the following list, in addition of providing an analytical expression 
for the nodal domain count, can also be implemented as a computer algorithm. We 
show that it performs as efficiently as the labeling algorithm. 

The counting methods that we present here are aimed for the discrete counting 
of both discrete and metric graphs. In what follows, we assume that a vector f is 
associated to the vertex set with entries fi. The nodal domains are defined with 
respect to f . 

4.1. Method I. - Counting nodal domains in terms of flips. We define 
a flip as a bond on the graph which connects vertices of opposite signs with respect 
to a vector f. The sign vector of f, denoted by f, is defined by fi = sign(/i). For 
the time being, it is assumed that f has no zero entries. The general situation will 
be discussed later. We denote the set of flips on the graph by ^(f): 

(14) T{i)^{{u,v)eB\^U<0} . 
The cardinality of J-{i) will be denoted by F(f ). 

Lemma 4.1. The number of flips of a sign vector f, can be expressed as 

(15) F(f)=F(f) = i(f,Lf). 

Proof. Using: 



(16) (f,if) = ^E (/"-/«) 



V^U 



, 2 _ ) 4' if /t> and /„ have opposite signs 

(1' ) [fv ^ fu) \ ^ ^ 

0, ii fv and /„ have the same sign 



D 



Using the number of flips, one can get an expression for the number of nodal 
domains: 

Theorem 4.2. Given a connected graph Q on V vertices, B bonds (and r 
cycles) and a vector f , then the number of nodal domains of f is: 

(18) v{^) = \C^,M) + V-B + m = J(f,Lf) - [r-m) + 1 

where /(f) is the number of independent cycles in Q of constant sign (with respect 
to f ) . The second equality above is based on equation (fTl . 

Proof. Let us remove all the flips from the graph. We are now left with a 
possibly disconnected graph Q. There is a bijective mapping between components 
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of Q and nodal domains of Q. Hence, the number of components in Q is equal to 
the nodal domains count of Q with respect to f . Let the number of nodal domains 
in Q be denoted by v{{). Using (nj), it is clear that for the i*'' component (where 
i — 1,2,..., J^(f)): ri = Bi — Vi + 1 where r.i, Bi and Vi are the number of cycles, 
bonds and vertices of the i*'' component, respectively. It is also clear that all the 
cycles in Q are of constant sign, since there are no flips in Q. Thus, by our notation 
ri = li. Let us sum over the components: 



(19) 



m = J2^. = ^(i?. - y, + 1) = (B - F(f)) ^V + z.(f ) 



1=1 



Combining (151 with (19 1, we get (18) 



D 



(18) is valid only for vectors f with no zero entries. In order to be able to 
handle a zero entry in f , we must perform a transformation on the graph. For 
a strong nodal count, we simply delete all the zero vertices along with the bonds 



connected to them from the graph, and then apply (18 1 on the new graph (with 



the new Laplacian). For a weak nodal count we replace each zero vertex by two 
vertices, one positive and one negative (not connected to each other), and connect 
them to all vertices which were connected to the original one. Now we can apply 



(18), and get the desired result. Notice that this correction fails in the case of a 



zero vertex whose neighbors are of the same sign (in this case an artificial nodal 
domain is added). However, the situation above can not occur for an eigenvector 
of a discrete Laplacian. This way of handling zero entries can be adapted for the 
following counting methods as well and will not be repeated in the sequel. 



Using (18), we can write some immediate consequences: 



(20) 
(21) 



F(f„) + 1 - r < ;/„ < F(f„) + 1 
n-r- 1 < F(f„) < n + r- 1 



(20 1 results from the obvious fact that < I < r, while (21) is a consequence 



of Courant's nodal domains theorem and Berkolaiko's theorem which states that 
n ~ r < Vr, ■ 



In order to make use of (18), one must compute l{{) which is not given explicitly 



in terms of f. Thus, it cannot be considered as an analytic counting method, nor 
does it offer computational advantage (There is no known eflicient algorithm which 
counts all the cycles of constant sign with respect to f). However, it offers a useful 
analytical tool for deriving other results, and it makes a useful connection between 
various quantities defined on the graph. 

4.2. Method //. Partition function approach. Foltin derived a parti- 
tion function approach to counting nodal domains of real functions in two dimen- 
sions [32 . It can be adapted for graphs in the following way: Each vertex, i, is 
assigned an auxiliary "spin" variable Si where Si = ±1 (a so called Ising-spin). 
Thus, given a certain function f on the graph, each vertex is assigned with two 
"spins" Si and fi. Let s denote the auxiliary spin vector: s = (si, S2, . . . , sv). 
Foltin introduced a weight to each configuration of the spins model. It assigns the 
value 1 to configurations in which all spins Si belonging to the same nodal domain 
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(with respect to f ) are parallel, while spins of different domains might have different 
values. The weight reads: 



^_i+M.i 



(22) w{{,s)= H 

It can be easily checked that this form satisfies the requirements stated above: The 
weight can take the values one or zero. It is one if and only if each factor in the 
product is equal to one. A certain factor is one, in either one of the two cases: if 
fi 7^ fj ihj belong to different domains) - this allows the Ising-spins in different 
domains to be independent of each other. The second case is if i,j are in the same 
domain, fi — fj and the corresponding Ising-spins are equal, Si — Sj. Let us now 
sum over all possible spin configurations {si} to get the partition function 

(23) Z(f)^^ti;(f,s). 

For the configurations whose weight has the value one, the spins have equal signs 
over each nodal domain and different domains are independent of each other. Hence, 
the total number of such configurations is: 

(24) Z(f) = 2'^(^) , 

where i>{{) is the number of nodal domains of the vector f. The nodal domains 
count is: 

(25) :/(f) = -'-lnZ(f) w 1.441nZ(f). 

The partition function approach provides an explicit formula for the number of 
nodal domains, and therefore it belongs to the analytic and not to the algorithmic 
counting methods. As a matter of fact, it is highly inefficient for practical computa- 
tions. It involves running over all possible spin configurations {sj}, where Si = ±1. 
There are 2^ different configurations, and as V increases the efficiency deteriorates 
rapidly. 

The partition function approach can be used as a basis for the derivation of 
some identities involving the graph and a vector f defined on it. It is convenient to 
introduce the following notations: 

1 + Mj 



fh 



o-fe 



2 

1 " S^Sj 



2 

Where f and s are as before, and h is an undirected bond. We generalize the 
partition function by introducing a new parameter x into the definitions 

(26) w{i,s;x) = ]J [1 - </?bcrfcx] 

beB 

(27) Z{i-x) = Y,w{i,s-x)=Y,\{{\-ipwbx) 

{s} {s} fcee 

(28) Z(f;l) = 2''(^) 
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where x can assume any real or complex value. At a; = 1, the generalized partition 



function is identical to ( 23 1 



Let us now perform the summation over all the vectors s, and compute the 
coefficient of x'^' . To get all the contributing terms we have to sum over all choices 
of k brackets from (27), in which x appears. Non vanishing contributions occur 
whenever both iph and at are equal to one. Since we are only summing over s, we 
only need to check when ct;, = 1. This happens if and only if the s vector has a 
flip on the bond b. Since we choose k brackets (which is equivalent to choosing 
k bonds), we need to count how many s vectors have flips on all these k bonds. 
The signs of those s vectors on bonds which are not contained in this choice of k 
bonds are irrelevant. If we observe the choice of k bonds (5i, . . . , 6^), we notice 
that each connected component, within this choice, contributes a factor of 2, since 
the symmetry of turning each plus to minus and vice versa, does not change the 
flip properties. Using (IT]) we see that the number of connected components with 
respect to the choice of k bonds is: Co{bi, . . . , bk) = V — k + r{bi, . . . , bk), where 
r(bi, . . . , bk) is the number of independent cycles that are contained in this choice. 
Finally we notice that a cycle of odd length cannot have flips on all of its bonds, 
so we will not sum over choices of 5i, . . . , 6^ which contain a cycle of odd length. 
Thus, the sum over s yields: 



fe=0 bi,...,bkeB 2=1 



(29) 



2"E 



fe=0 



til 
2^ 



6i,...,fcfcee i=i 



Where 



stands for summation on all the possibilities to choose k bonds 



■^^ 6i....,6fceB 

bi,. . . ,bk G B such that the subgraph they form do not contain an odd cycle. We 
can now derive some immediate properties of the generalized partition function. 
To start, we compute the leading four derivatives at x = to demonstrate the 
counting techniques involved. Some more effort is required to compute the higher 
derivatives. 



(30) Z(f; 0) 

(31) Z(i)(f;0) 

(32) Z(2)(f;0) 

(33) Z(3)(f;0) 



-2^-iEVb = -2^~'(S-F(f)) 
fcee 

2!2^-E' fl(^^.) = 2"-^P7^'^ 



bi,b2&B i=l 



-3! 2^-3 E n(^0 = -3!2 

bi,...,b3eB i=l 



V-3 



B-F{{) 
3 



C. 
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(34) 



zW(f;0) = 4! 2^-4 Y! W^"^) 

bi,...,bieB i=l 

r/B-F(f) 

^"il 4 
B-F(f) 
4 



4! 2^-4 <j 2G 
4! 2^-4 ' 



C4 — Cq 



+ C4 — C3 



Where C3 is the number of triangles of constant sign, C4 is the number of cycles 
of length 4 of constant sign, and C3 — 6*3(5 — F(f) — 3) is the number of choices 
of 4 non-flips bonds which contain a triangle. In the evaluation of the function and 
its flrst three derivatives we have used the identity Vn < 3 r(6i, . . . , 6„) — when 
&i, . . . , 6„ contain no odd cycles. 

The partition function is a polynomial of degree less than or equal to the 
number of bonds. For a graph Q with B bonds and Co connected components, 
the polynomial is of degree B if and only if Q is bipartite and the function f is of 
constant sign on each connected component. This follows from two observations: 
first, a well known theorem in graph theory states that a graph Q is bipartite if 
and only if it contains no cycles of odd length. Thus we can choose all the bonds 
in Q without having an odd cycle contained in this choice. Second, unless f is of 
constant sign on each connected component, we will encounter a flip and hence the 
multiplication over all ipt- will vanish. If Q is bipartite, the coefficient of x^ is 2*-^°. 

While the value of the polynomial at x = 1 has an immediate application 



through (28 1, we can evaluate it for other values. Let us choose f to be a vector of 
constant sign. This way (fti = 1 for all i = 1, 2, ... , B. Hence (29 1 reduces to: 



(35) 



2"E^?I E' ^''''- 

fc=0 \6i,...,6fcG8 



b.) 1 ^k_ 



On the other hand, (27 1 equals: 
(36) 



J2llil-a,x). 

{s} beB 



Choosing x ~ —n+ 1 where n G Z and using the fact that (35 1 and (36 1 are equal, 
we get: 



(37) 



;vE' 



2^ 



,F(s) 



{s} 




bi,....bkeB 



For random vectors s, uniformly distributed, the left hand side can be interpreted 
as the average over this ensemble of the quantity: n^'^^K Let us choose two special 
values for n. If we choose n = —1, the left hand side is just the difference between 
the probability that a random vector s on the graph has an even number of flips, 
and an odd number of flips. The right hand side is: 



(38) 



E(-i)' E 2"^^'^ 

fc=0 bi,...,bk€B 



-M) 



NODAL DOMAINS ON GRAPHS - HOW TO COUNT THEM AND WHY? 



13 



For n = 2, we get: 

(39) ^E2^^^^=E4 E' 2^ 



{«} 



2fc 
k=o fci,...,fcfcee 



{bi,....M 



An example of the use of the polynomial could be in proving that on a tree, the 
probability of an even number of flips is equal to the probability of an odd number 



of flips. We use Eq. (38 1 and see that this difference of probabilities is equal to: 

Other possible identities we can derive are for the complete graph, Ky- A function f 
of constant sign induces one nodal domain on Ky , while any other function induces 
two nodal domains. Therefore, for a; = 1: 



(40) 
(41) 



2"E 



(-1)^ 
ofe 



Y^ 2'-(^i 



M 



■iV 



B 



E E 



(-1)'^ 



2^-2 ^^ -^ 2^ 
fe{±i}''' fe=o 



fe=0 bi,...MkeB 

E' fn^^.(f))2^^^^' 

6i,...,6fcee \i=i / 



,bk) 



Where fbii^) — 2 ^'^^ ^i ~ {''^j'")- Equivalently, running over all the functions 
f (including those of constant sign) we get: 



(42) 2^ E E 



(-1)* 



fe{±i}^fc=o 



E' n^^.(f) 2^^^^' 



6i,...,6fcee \i=i 



= 4(2^ - 1) 



4.3. Method III. — Breaking up the graph. We begin again by deleting 
all the flips from the graph Q. This way we are left with a (possibly) disconnected 
graph, Q in which each connected component corresponds uniquely to a nodal 
domain in the original graph. The connectivity matrix C and the discrete Laplacian 
L oi Q are given by 



(43) 
(44) 



Gi 



C\j ■ 



— L/,;. 



1 + hfj 



Oij / ^ Cik 



fc=l 



Where f is the sign vector, and it is assumed for the moment that none of the 
entries of f vanish. 

We now make use of the theorem which states that the lowest eigenvalue of the 
Laplacian is with multiplicity which equals the number of connected components 
in the graph. Therefore, finding the nodal domains count reduces to finding the 
multiplicity of zero as an eigenvalue of L. An analytic counting formula can be 
derived by constructing the characteristic polynomial of L: 



(45) 



det(A/y - L) 
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The multiplicity of its eigenvalue provides the nodal domains count: 
(46) iy(f) = lim A— lndet(A/v - L) . 



This method of counting, which provides the analytical expression (46 1 for the 
nodal count, is also the basis for a computational algorithm which turns out to be 
very efficient. It relies on the efficiency of state of the art algorithms to compute 
the spectrum (including multiplicity) of sparse, real and symmetric matrices. To 
estimate the dependence of the efficiency on the dimension V of the graph, we have 
to consider the costs of the various steps in the computation. The construction of 
the matrix L, takes 0{V^) operations, and storing the information requires 0{V^) 
memory cells as well. It takes 0{V°') operations to find all its eigenvalues where 
a ~ 2.3 (and at worst case a ~ 3) [29j . In figure fl] this polynomial dependence is 
shown for graphs of two different connectivity densities, with y equals 0.5 and 5. 
In this figure the logarithm of the time needed to find all eigenvalues of L (defined 
for a random vector), is plotted against the logarithm of the number of vertices. 
The slope which is the exponent of the polynomial dependence is smaller than 3. 
The eigenvalues in these two examples were attained using the Matlab command 
eig. As will be shown below, there are more efficient ways of finding the spectrum 
of sparse, real and symmetric matrices. Thus, the efficiency stated above can be 
improved for graphs with sparse Laplacians. 

Finally, we would like to show that the present method can be applied for 
counting nodal domains of functions defined on two dimensional grids and that 
its efficiency is comparable to that of the commonly used HK algorithm. Given 
a function f on a two dimensional domain, we have to compute its values on a 
rectangular grid with ^/Vx \/V points. The HK algorithm counts the nodal domains 
in 0(y) operations |30j . Using our method, we consider the rectangular grid as a 
graph with V vertices. Assuming for simplicity periodic boundary conditions, the 
valency of all the vertices is 4. The corresponding L matrix is a y x ]/ matrix 
which is sparse (as long as F ^ 4), and due to the periodic boundary conditions it 
takes the explicit form: 

(47) iij- = 4(5ij- — (5ij-i — ^ij+i — ^i-v,2 — ^i+v,j ■ 

Thus, storing L takes only 0(y) memory cells, and constructing L takes 0(y) 
operations. We mentioned above that for a general real symmetric matrix, the 
number of operations needed is 0(y"), where a ~ 2.3 and at worst case a ~ 3. 
However, the sparse nature of L significantly simplifies the problem. The most 
well-known eigenvalue method for sparse-real-symmetric matrices is the Lanczos 
method. In addition, in recent years, new efficient algorithms were discovered for 
the same problem. In |31j . it is proven that finding the eigenvalues of a sparse 
symmetric matrix takes only 0(y) operations. Combining the costs, we find that 
it takes our algorithm 0{y^ operations in order to compute the nodal domains 
count, and therefore it is comparable in efficiency to the HK algorithm. 
As mentioned earlier, the labeling algorithms also display linear efficiency f[46 .j47j). 
The labeling algorithms have the advantage that they are simpler in a sense, and 
that they are implemented quite easily as computer programs. In addition, the la- 
beling algorithms maintain their linear efficiency even for graphs with dense Lapla- 
cians. It is worth mentioning, however, that our algorithm has the advantage that 



it provides an analytic expression of the nodal domains count ( 46 1 
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Figure 1. The time it takes to compute the spectrum of L as a 
function of the number of vertices, for two different connectivity 
densities: t7 = 0.5 in the upper figure and ^ = 5 in the lower 
figure. 



4.4. Method IV. — A geometric point of view. The counting method 
proposed here uses a geometric point of view which starts by considering the V 
dimensional Euclidean space, and dividing it into 2^ sectors using the following 



^ ' f->^ ' /-,^ ' 1 TTT K J-ll-J-l fl^ 



el 



construction. Consider the 2^ vectors e'"^ = (e]^" , 63 , . • . ,ey ) where 
{1,-1}, a = 1,2,..., 2^. A vector xe M^ is in the sector a if x-e'") = X^Li 
In two dimensions, the sectors are the standard quadrants. We shall refer to the 
vectors e^") as the indicators. 



Xi 



16 RAM band", IDAN OREN" AND UZY SMILANSKY"'^ 

Given a graph Q with V vertices, we partition the 2 indicators into disjoint 
sets: 7„ = {e^") : (^(e^"^) = n} where iy{e'^°'^) denotes the nodal domains count of 
the indicator e^"^ with respect to Q. As shown before max{n |7n7^0} <V — x + 2, 
where x is the chromatic number of Q, and also some of the 7^5 might be empty. 

Let f be a vector with non-zero entries defined on the vertex set of Q. Then, 
the main observation is that i^(f) = n if and only if: 

(48) Y. ^((e",f)-E|f. 1) = 1, 

where, 

/A,^\ i/ \ T e . 2; f 1, if a; = 

(49 5{x) = hm - sm - = <^ ' . 

e^o X e [ 0, if a; ^t 

and (e", f) is the dot product of e" and f. In other words, by finding the sector to 
which f belongs and knowing from a preliminary computation the number of nodal 
domains in each sector, one obtains the desired nodal count. Thus, the present 
method requires a preliminary computation in which the sectors are partitioned 
into equi-nodal sets 7„. This should be carried out once for any graph. Therefore 
the method is useful when the nodal counts of many vectors is required. In several 
applications, one is given a vector field (of unit norm for simplicity) f £ §^^^ which 
is distributed on the {V — l)-sphere with a given probability distribution p{i), and 
one is asked to compute the distribution of nodal counts, 

(50) P{n)= I p(f)(5(i^(f)-n) d^-if . 

In such cases, the preliminary task of computing the equi-nodal pays off, and one 
obtains the following analytic expression for the distribution of the nodal counts. 

P{n) = / p(f) d^-if ^ 6 [ (e", f) - f] |f,| 1 = 
f p(f)5(l-|fp)d^f ^ 5((e",f)-5]|f,|l 



(51) 



where: p(f) = p(^) • 2|f|. 



.|f|- 



(51 1 can also be formulated as: 



(52) 



p{n) = Y.I P(^)^(i-ir)d^f E ^f(-"'f')-E/«l 
E /%(f'')<5(i-ifnd^f E ^f(e",f")-E/.) 

ef^etr, e°e7„ \ 1=1 / 



Where / — /f.>Q means integration on the first sector (the vectors with all entries 
positive) and f = (/le^', /2e^, . . . , fvey)- 



Equations (51 1 and (52) are the general equations governing the nodal domains 



count distribution. In order to make further progress, we need to specify the dis- 



tribution from which f is taken. This means that we need to specify p{f) in (51) 
for example. Let us discuss two examples: 
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A uniform dist ribu tion over the V — I dimensional sphere: In this case, we can 



solve Equation (51) and get that P{n) ~ '-^. Note that for a tree, we can solve 



this problem by other means. Using (18 1, we see that for a tree, the number of 



nodal domains is equal to the number of flips plus one. Since f is taken from the 
uniform distribution, then the probability of a flip is half. The number of flips in 
a vector f is thus a binomial variable: F(f) ~ Binomial{N,p) with N = V — 1 is 
the number of bonds, and p = |. For large enough V this approaches the Gaussian 
distribution: F(f) ^ Gaussian{^,a^) with /i = ^^^ and cr^ ~ ^^^t^- From this 
result we can infer that: 

(53) P{n) « , exp (^^^^ ^^) 

For the other extreme, the complete graph, Ky, the only possible nodal domains 
counts are one and two |28j . The vectors which yield a nodal domains count of one 
are vectors of constant sign. All other vectors yield a nodal domains count of two. 



Indeed, using (51) or (52) it is easy to be convinced that for the complete graph, 
7i = 2 while 72 = 2^ — 2. All other 7„'s are empty. 

Micro- canonical ensemble: In this case the vectors f are uniformly distributed on 
the energy shell, where we can also deflne a measurement tolerance factor, A: 

,,0 ,n _ Mg-|(f,^f)-A|)^(l-|fp) 

In order to make use of this ensemble, further work must be done, for example, a 
natural way to order the functions of the ensemble. 



5. The resolution of isospectrality 

There are several known methods to construct isospectral yet different graphs. 
A review of this problem for discrete graphs can be found in |33j . The conditions 
under which the spectral inversion of quantum graphs is unique were studied pre- 
viously. In |38L I39j it was shown that in general, the spectrum does not determine 
uniquely the length of the bonds and their connectivity. However, it was shown in 
^35j that quantum graphs whose bond lengths are rationally independent "can be 
heard" - that is - their spectra determine uniquely their connectivity matrices and 
their bond lengths. This fact follows from the existence of an exact trace formula 
for quantum graphs |40], 141] . Thus, isospectral pairs of non congruent graphs, must 
have rationally dependent bond lengths. An example of a pair of metrically distinct 
graphs which share the same spectrum was already discussed in |35j . 

The main method of construction of isospectral pairs is due to Sunada |34j . 
This method enabled the construction of the flrst pair of planar isospectral domains 
in R^ |36j which gave a negative answer to Kac's original question: 'Can one hear 
the shape of a drum?' |37j . Later, it was shown that all the known isospectral 
domains in M? |42L I43j which were also constructed using the Sunada method have 
corresponding isospectral pairs of quantum graphs |44j . An example of this corre- 
spondence is shown in flgure[2] As mentioned in the introduction, it is conjectured 
[3 I19j that nodal domains sequences resolve between isospectral domains. For flat 
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Figure 2. (a) Planar isospectral domains of the Ta type, (b) Re- 
ducing the building block to a 3-star. (c) The resulting isospectral 
quantum graphs. 



tori in 4-d, this was proven [8j. We present here three additional known results for 
the validity of the conjecture for graphs. 

The first result is for the quantum graphs shown in figure p[c) . Both graphs of 
this isospectral pair are tree graphs and therefore have the same metric nodal count 
fJ-n = n |15j . This demonstrates the need to use the discrete nodal count in order to 
resolve isospectrality in this case. Indeed numerical examination of this case shows 
that for the first 6600 eigenfunctions there is a different discrete nodal count for ~ 
19 % of the eigenfunctions. Similar numerical results exist for two other pairs of 
isospectral graphs that are constructed from the isospectral domains in |42L I43j . 
The exact results are described in 19 . 

Another result is also in the field of quantum graphs |19j . The graphs in figure 
[3] are the simplest isospectral pair of quantum graphs known so far. The simplicity 
of these graphs enables the comparison between the nodal counts of both graphs. 
It was proved that the nodal count is different between these graphs for half of 
the spectrum. This result was proved separately for the discrete count and for the 
metric count. The proof does not contain an explicit formula for the nodal count 
but rather deals with the difference of the nodal count between the graphs averaged 
over the whole spectrum. 

Examining the nodal sequences for the graph II for various values of the length 
parameters a, b, c, we observed that the formula 

(56) M^^ = --2-2(-l)'^^"' 
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I II 

Figure 3. The isospectral pair with boundary conditions. D 
stands for Dirichlct and N for Neumann. The bonds' lengths are 
determined by the parameters a,b,c 

reproduces the entire data set without any flaw H . Assuming it is correct (which 
is not yet proved rigorously) we first see that it provides an easy explanation for 
the previously discussed result regarding the resolution of isospectrality for this 
pair. For rationally independent values for the parameters a,b,c one gets that 



^J■ 



II 



y^ n for half of the spectrum. Combining this with /if^ — n (since graph / is 
a tree) we see again that for half of the spectrum the nodal domain sequences are 



different. Expression (56 1 is a periodic function of n with period proportional to 
the length of the only loop orbit on the graph (the length is measured in units of the 
graph's total length). It can be expanded and brought to a form which is similar 
in structure to a trace formula where the length of this orbit and its repetitions 
are the oscillation frequencies. A similar trace formula for the nodal counts of the 
Laplacian eigenf unctions on surfaces of revolution was recently derived [6] . 

Finally, we direct our attention to discrete Laplacians. It was recently shown 
|28j that if Q and H are two isospectral graphs where one of them is bipartite 
and the other one is not, then their nodal domains count will differ. Without loss 
of generality, let Q he a, bipartite graph and Ti. a non- bipartite one, then for the 
eigenvector of the largest eigenvalue, the nodal domains count are different: for 
G, Vv — ^; while for Ti, v^ < V. The proof of this theorem is based on another 
interesting result derived in [28j: Denote by fy the eigenvector corresponding to the 
largest eigenvalue of the Laplacian of a connected graph Q. Then i'{iv) = vg = V , 
if and only if Q is bipartite. Figure 4. illustrates this result. 



6. Summary and open questions 

In spite of the progress achieved recently in the study of nodal domains on 
graphs, there are several outstanding open problems which call for further study. 
We list here a few examples. 

Of fundamental importance is to find out whether there exists a "trace formula" 
for the nodal count sequence of graphs, similar to the one derived in [6j for surfaces 
of revolution. The closest we reached this goal is for the graph II in the previous 



section, where ( 56 1 could be expanded in a Fourier series. However ( 56 1 was deduced 



numerically but not proved. Once a nodal trace formula is available, it could be 



This result was obtained with A. Aronovitch. 
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Figure 4. The upper figure presents a pair of isospectral graphs 
taken from [26 . Graph G, on the right is bipartite, whereas graph 
H , on the left, is not. The lower figure presents the nodal domains 
count, v{^n) vs. the index n. 



compared to the spectral trace formula |41j and might show the way to prove or 
negate the conjecture that counting nodal domains resolves isospectrality O I19j . 

The conjecture mentioned above can be addressed from a different angle. One 
may study the various systematic ways to construct isospectral pairs and investigate 
the relations between the construction method and the nodal count sequence of the 
resulting graphs. Such an approach worked successfully for the isospectral graphs 
presented in figure |3] [19' . 

Another open question which naturally arises in the present context: Can one 
find graphs whose Laplacians have different spectra but the nodal count sequences 
are the same? A positive answer is provided for tree graphs |15j . Are there other 
less trivial examples of "isonodal" yet not isospectral domains? 

It follows from Berkolaiko's theorem |27j that the number of nodal domains 
(both metric and discrete) of the n*'' eigenfunction is bounded in the interval [n — 
r, n] . We can thus investigate the probability to have a nodal count Vn = n — f (for 
< f < r). This probability, which is defined with respect to a given ensemble of 
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graphs, is denoted by P{r). It is defined for discrete graph Laplacians as: 

(57) Pir) = ^{if{l<n<V: ly^ ^ n ^ r}) . 

The corresponding quantity for metric Laplacians is: 
N{K) = {# {n : kr, < K}) 
Pir;K) ^ _i_^ (#{„<7V(X) : ^„ = „_f}) 

(58) P{r) = lim P{P;K) 

AT— >oo 

Here, ( ) stands for the expectation with respect to the ensemble. New questions 
arise from the investigation of the relation between the connectivity of the graph 
and the nodal distribution P(f) . Can one use the information stored in P{f) to gain 
information on the graphs e.g., the mean and the variance of the valency (degree) 
distribution of the vertices in the graphs? 

Many of the results we have presented, have analogues in Riemannian manifolds 
(which in most cases, were discovered earlier) - for example, Courant's theorem 
was originally formulated for manifolds. One can search for other analogues, and a 
good example is the Courant-Herrmann Conjecture (CHC). For manifolds the CHC 
states that any linear combination of the first n eigenfunctions divide the domain, 
by means of its nodes, into no more than n nodal domains. Gladwell and Zhu |45j 
have shown that in general there is no discrete counterpart to the CHC. However, 
we can still ask for which classes of graphs does the CHC hold? 
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